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Abstract 



This Master of Science thesis treats the hypercharge axion, which is a hy- 
pothetical pseudo-scalar particle with electroweak interactions. 

First, the theoretical context and the motivations for this study are 
discussed. In short, the hypercharge axion is introduced to explain the 
dominance of matter over antimatter in the universe and the existence of 
large-scale magnetic fields. 

Second, the phenomenological properties are analyzed and the distin- 
guishing marks are underlined. These are basically the products of photons 
and Z'^s with high transverse momenta and invariant mass equal to that of 
the axion. 

Third, the simulation is carried out with two photons producing the 
axion which decays into Z^s and/or photons. The event simulation is run 
through the simulator ATLFAST of ATLAS (A Toroidal Large Hadron Col- 
lider Apparatus) at CERN. 

Finally, the characteristics of the axion decay are analyzed and the crite- 
ria for detection are presented. A study of the background is also included. 
The result is that for certain values of the axion mass and the mass scale 
(both in the order of a TeV), the hypercharge axion could be detected in 
ATLAS. 
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Introduction 



Ever since the discovery of antimatter it has been a mystery that almost all of 
our experiments in particle physics are symmetric in matter and antimatter, 
but yet the universe seems to be constituted entirely of matter. Some general 
conditions were outlined by Sakharov in 1967, but the problem itself remains 
unsolved. 

The first possible explanations were based on the grand unified theories 
(GUTs) in which the asymmetry was generated very close to big bang. These 
theories allow baryon-to-lepton decay, which could generate an asymmetry. 
One drawback is that the theories cannot be tested without tremendous 
amounts of energy (> 10^^ GeV), far beyond our reach. 

Later on, theories evolved that could explain the generation of the asym- 
metry without GUT theories. Most of them suppose that the asymmetry 
was created around the electroweak phase transition at T ~ 100 GeV when 
the electroweak symmetry was broken. Supersymmetric theories can offer 
possible explanations under certain conditions, but the standard model itself 
seems incapable to produce the required asymmetry. 

A theory proposed by Brustein and Oaknin is particularly appealing for 
several reasons. It suggests simply that the introduction of a scalar field 
could create the required asymmetry through coupling to the hypercharge. 
As a bonus, the theory could also explain the existence of large magnetic 
fields in cosmic plasmas. The pseudoscalar was named hypercharge axion 
from its couplings to hypercharge. 

In summary, the hypercharge axion can possibly explain the dominance 
of matter over antimatter in the universe. In other words, why we exist! 

Finally, the hypercharge axion is expected to have a mass in the TeV 
range, making it possible to detect in the ATLAS detector at the large 
hadron collider (LHC) which is now constructed at CERN, Geneva, and 
will be operational by 2005. This is the subject of this masters thesis. 

Some references on particle physics are Nash [||, Peskin and Schroeder 
P], and the CERN photo-gallery: http: / / press. web. cern.ch /Press /Photos / . 
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Chapter 1 

Background and Motivations 



This chapter treats the background to my research, why it may be of interest 
to study the axion and what has been done so far. Briefly, the axion could 
be the reason for the mysterious domination of matter over antimatter in the 
universe. Before introducing the axion some theory is presented to situate 
it in its context and explain how the axion can be the reason behind the 
domination of matter. 

This chapter is divided into three major parts. The first one addresses 
some general questions, such as how we can be sure that the universe really 
has more matter than antimatter, and some very general conditions that 
must be fulfilled for the asymmetry to exist. These conditions are called the 
Sakharov criteria and they must be fulfilled for statistical reasons. 

The second part of this chapter is devoted to baryogenesis, i.e., the pro- 
cess by which the asymmetry is created. The subject is very complicated, 
so I will content myself with a rather qualitative description of the differ- 
ent phenomena. Among these are the electroweak phase transition, which 
took place about lO"^'^ seconds after the big bang; the role of hypermag- 
netic fields in the phase transition; sphalerons which counter asymmetries 
and the Chern-Simons number, which is directly proportional to the baryon 
asymmetry. 

The third part treats the axion itself, with the properties proposed by 
Brustein and Oaknin A description of how the axion could amplify 
the hypermagnetic fields is included as well as the couplings and branching 
ratios. A discussion of where the hypercharge axion could appear is also 
provided. Its properties are compared to those of particles in the most 
popular extensions of the standard model: the minimal supersymmetric 
model and superstring models. 
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1.1 Matter- Antimatter Asymmetry 



1.1.1 How do we know that the Universe is Asymmetric? 

While our theories in particle physics are very nearly symmetric in matter- 
antimatter (almost all relevant experiments in particle physics show this 
behavior), the universe itself is far from being so. We have several reasons to 
believe that the universe is mainly constituted of matter. The most obvious 
reason is that we have explored our solar system without annihilating, which 
would not have been the case if there had been antimatter in significant 
abundance in our solar system. 

To conclude that the universe at large is made of matter is not as ob- 
vious. What if another part of our galaxy, or another galaxy nearby were 
constituted of antimatter? This possibility is improbable since cosmic radi- 
ation in that case should contain antimatter. Cosmic radiation does contain 
some antiprotons, but that fraction (10~^) is well explained by the reaction 
p+p ^ 3p+p, which has nothing to do with the eventual existence of larger 
chunks of antimatter. This is evidence that there is no large antimatter 
cluster in our relative proximity. 

As for more distant galaxies, we should see gamma-ray emissions from 
the interface between matter and antimatter galaxies. The gamma rays 
should have an energy up to that of a particle-antiparticle pair. These 
gamma rays should produce a detectable background, which is not observed. 
Causality is not an argument against this either, as we can see (almost) all 
parts of the universe, only at different times. The separation should have 
occurred rather early for any significant amount of antimatter to remain, so 
that we should observe gamma rays from these parts too. Thus, there is a 
negligible amount of antimatter on the scale of clusters. For the interested 
reader, references like Steigman Stecker and Cohen et al. Q are 
recommended. 

1.1.2 How did the Universe become Asymmetric? 

The first possiblility is that the universe started asymmetric. This is not a 
very interesting possibility because of its sterility. It does not suggest any 
new physics, nor any "real" explanations or predictions. Besides, it is in 
contradiction with inflation theory, which says that any initial abundances 
are diluted. 

The second possibility is that the universe started symmetric (or at least 
became symmetric at a rather early stage, e.g. during inflation), but the sym- 
metry was broken through some mechanism. For the baryon asymmetry this 
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mechanism is called baryogenesis. Note that even though baryon number 
and lepton number separately are violated, the difference B — L \s thought 
to be constant. 

The imbalance in the amount of matter and antimatter can basically be 
measured as the imbalance in the amount of baryons and antibaryons, since 
baryons (neutrons and protons) are the heaviest of the basic constituents of 
normal matter. Also, we cannot measure the lepton asymmetry because we 
know neither how many neutrinos there are in the universe, nor their mass. 
The asymmetry is often characterized by the baryon-to-entropy ( [pO| , page 
133) ratio: 

= ^ = « 10-1°, (1.1) 

s s 

which can be seen as the the ratio between the number of baryons and 
the number of photons (which come from the annihilation of matter and 
antimatter). 

In 1967, Sakharov |^ stated that there are three criteria that must be 
satisfied to explain the present matter-antimatter asymmetry: 

1. Non-conservation of baryons 

If the baryon number is conserved in all reactions then the observed 
asymmetry can only reflect asymmetric initial conditions. Hence, there 
must have existed some process that violated the baryon number. 
Grand unified theories (GUTs) have leptons and quarks in the same 
representation, and therefore allow decay of baryons to leptons. A pro- 
cess that has been searched for is proton decay to positron, p — > e^Ve, 
but no evidence for such processes has been found so far. However, 
baryon number violation can also be obtained from anomalies, which 
can occur without any extension from the standard model, see below. 

2. Charge (C) and Charge-Parity (CP) violation 

CP symmetry implies that if there are two identical decay channels, 
except that one involves matter and the other involves antimatter, 
then both happen with equal probability, e.g., P{e~p'^ — > niy) = 
P{e^p nv). In other words, if CP symmetry is violated and if 
you could contact someone in another part of the universe you could 
instruct him how to make an experiment to tell whether he is made of 
matter or antimatter. Without a preference for matter or antimatter, 
baryon non-conservation would work as much one way as the other, 
and matter and antimatter would have annihilated. 
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3. Thermodynamic nonequilibrium 

The universe must have evolved from a state of thermodynamic equi- 
hbrium to thermodynamic nonequilibrium. This is necessary because 
in thermodynamic equilibrium the net number density of baryons is 
equal to that of antibaryons. This can be illustrated as follows. In 
thermodynamic equilibrium the number density n{E) of a particle is 
given by 

n{E) oc e-^^^/^, (1.2) 

where ks is Boltzmann constant, T is the temperature and H is the 
Hamiltonian of the system. Now, due to the Charge-Parity-Time 



(CPT) theorem |1C] the masses of particles and antiparticles are the 
sameQ, which means that the Hamiltonians are the same. In other 
words, equation ( |1.2[ ) implies that in thermodynamic equilibrium, the 
number (density) of baryons is the same as the number (density) of 
antibaryons. You can also explain this condition with the fact that 
a difference in the number of particles and antiparticles decreases the 
randomness of the system, which is prohibited in thermodynamic equi- 
librium, as the second law of thermodynamics states that the entropy 
of a system must increase, or at best stay the same. 

One way of bringing the universe out of thermodynamic equilibrium 
is by a phase transition. However, the thermodynamic nonequilib- 
rium needs to be very strong to allow the current asymmetry between 
matter and antimatter. A second-order phase transition is charater- 
ized by continuously varying parameters, giving a weak thermody- 
namic nonequilibrium. A first-order phase transition has discontinu- 
ously varying parameters. To create the currently observed matter- 
antimatter asymmetry, a strong first-order phase transition is needed. 

In summary, the entire production of the asymmetry should go trough 
the following steps (assuming the simplest procedure): 

1. The universe started symmetric in matter and antimatter. 

2. Processes violating baryon number, C and CP symmetry came into 
action. 

3. The universe was (temporarily) brought out of thermodynamic equi- 
librium, allowing for an asymmetry between matter and antimatter to 
develop. 



^ The validity of the CPT theorem is very fundamental, but is nevertheless questioned, 
see for example Eades and Hartman |llt|. 
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4. Very shortly after (or during) the return to thermodynamic equiUb- 
rium the baryon number violating processes stopped (otherwise the 
return to thermodynamic equilibrium would once again have brought 
the universe to symmetry between matter and antimatter). 

5. The remaining matter and antimatter annihilated (into photons) until 
only the surplus matter remained. 

1.2 Electroweak Model of Baryogenesis 

In the minimal standard model, two of the Sakharov criteria are problematic. 
First, the high mass of the Higgs boson seems to preclude a first-order phase 
transition Q, which would have been a natural way to bring the universe 
out of thermodynamic equilibrium. Second, the CP violation is too weak 
to produce the observed baryon-to-photon ratio lO"^'^. This means that 
extensions of the standard model need to be considered, but before doing 
that an explanation of how the asymmetry can be created is appropriate. 

Baryogenesis should take place sometimes very early in the history of the 
universe, in the first fractions of a second, before the formation of compounds 
of quarks like baryons and mesons (note that "free" quarks also have a 
baryon number). 

As the baryogenesis must take place in thermodynamic nonequilibrium, 
there are two basic possibilities to investigate. First of all, the pure expan- 
sion of the universe provokes a thermodynamic nonequilibrium. However, 
at the time around the electroweak phase transition this has been estimated 
to be by far too weak to stop the symmetry from being restored. Second, 
the thermodynamic nonequilibrium can be caused by a phase transition, like 
the electroweak phase transition 10"^'^ seconds after the Big Bang. Thus 
the thermodynamic nonequilibrium must either be caused by a phase transi- 
tion, or occur at a very early stage at the scale of the grand unified theories 
(GUTs). In this study, the electroweak phase transition will be treated in 
further detail. 

Some aspects of the GUT thermodynamic nonequilibrium can be found 
in Kolb and Wolfram |^]. 

1.2.1 Analogy with a Pendulum 

For several of the concepts introduced in the following sections, an analogy 
with a mechanical pendulum can be of great help. The Lagrangian of the 
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pendulum is: 



L = 



\mf'6'^ - mgl{l — cosO) 



(1.3) 



where m is the mass and I is the length of the pendulum (6 is the angle of 
inclination of the pendulum with respect to the vertical and g « 9.8 m/s^ 
is the standard acceleration of gravity). 

The system is periodic and the periodicity can be labeled by an integer 

n: 



Below is a list of concepts that have anlogues for the electroweak phase 
transition: 

Winding number The number of rotations, n. 

Vacuum The solution where the pendulum is at rest = 0). 

Sphaleron The solution where the pendulum is at its highest point = vr). 
Sphaleron is also often used as a label of a rotation, passing the 
sphaleron point, going to an adjacent winding number n. The energy 
of the sphaleron point is mgl{l — cosvr) = 2mgl. 

Instanton The tunneling of the pendulum from one winding number to 
the adjacent one. (This is a quantum mechanical effect and does not 
really make sense in this mechanical context, but it is included for 
completeness.) 

Thermal energy A set of pendula having a randomly distributed energy, 
but with the mean energy equal to the thermal energy. 

If the thermal energy becomes comparable to the sphaleron energy, it be- 
comes possible for thermal transitions over the energy barrier to occur. The 
rate of these transitions can be calculated to be 



so that when T ^ Egph the transitions between different vacua (different 
winding numbers) happen unhindered. 

Another analogy between the pendulum and quantum field theory is that 
for low T, i.e., ^ w 0, we can approximate cos9 1 — 9^/2, which is the 
same thing as perturbation theory, which also is a kind of Taylor expansion 
around a minimum in the potential. However, in both cases information 
about the periodic structure of the system is lost. 



On = 2mr. 



(1.4) 



r(r) oc e-^'p'^/'^ 



(1.5) 
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1.2.2 Thermodynamic Nonequilibrium 

This is the era when the Higgs first acquired a vacuum expectation value 
(vev), so the and acquired masses. The temperature of the universe 
was at this stage of the order of 100 GeV. 

As mentioned earher, the phase transition has to be of first-order to 
induce a thermodynamic nonequihbrium strong enough. In fact, it is not 
even sufficient with a first order phase transition, but a strong first order 
phase transition is required in the simplest model. 

The dynamics of the phase transition is obviously very important, and 
the imagined procedure is in analogy with that of pressurized boiling water. 
Close to a critical temperature Tc, "bubbles" of the other phase starts to 
form. These are regions of space where the electroweak symmetry is broken. 
For a start these bubbles die out as fast as they appear, but then they grow 
fast enough to form nucleates (called critical bubbles) and they expand and 
eventually fill the entire universe. At the bubble walls there is a significant 
departure from thermodynamic equilibrium. 



1.2.3 Nonconservation of Baryon Number 

Classically^ the Lagrangian conserves both the baryon number and the lep- 



ton number. However, there are triangle anomalies (see figure 1^), which 



add radiative corrections to the Lagrangian. The theory can be renormal- 
ized with conservation of either the vector or the axial current, but not 
both. We know very well that the vector current is conserved (by conserva- 
tion of charge), which means that the axial current is not conserved. Now 
the integral over time of the divergence of the axial current (along with the 
conservation of the vector current) implies a change in baryon number 

However, the radiative corrections from the anomalies are non-perturbative 
and at the temperature of the universe today, the corrections are extremely 
small. In the case of the pendulum above, they correspond to a thermal ex- 
citation such that the pendulum would pass its highest point, the sphaleron, 
or they correspond to instanton tunneling. 

As in the case of the pendulum above, there is also a winding number 
associated with the Lagrangian, the Chern-Simons winding number, Nqs- 
The relation between the Chern-Simons number and the change in baryon 



^ Classically means with only perturbative calculations. 

^ Compare with classical electron current V • J = AQ, where J is the current, t is 



time and AQ the change in charge. 
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Figure 1.1: The triangle anomaly contributing to the baryon and lepton 
number currents, and thereby the violation of baryon and lepton number. 



number can be expressed as: 

= ANcs = nflNcsih) - Ncsito)], (1.6) 

where nf is the number of families and Ncsit) is the Chern-Simons (wind- 
ing) number at time t. Hence, when the Chern-Simons number is changed 
by one, nine quarks, (three color states for each generation) and three lep- 
tons are created, giving AB = 3 and AL = 3. This change in Chern-Simons 
number is vehicled by either a sphaleron process or an instanton. Without 
pretending to know the details, I give the basic procedure as follows. 



The electroweak potential can be sketched as in figure |1.2| . There are 
several different vacua, like in the case of the pendulum, and the different 
vacua have different Chern-Simons number (like n in the case of the pendu- 
lum). A transition between different vacua (and therefore different Chern- 
Simons number and hence a change in baryon number) can be obtained in 
two different ways. First, a tunneling event may occur, which is called an 
instanton. The probability of tunneling depends on the thermal energy, i.e., 
how "thick" the potential barrier between different Chern-Simons numbers 
is. Second, if the (thermal) energy is high enough, it is possible to pass 
above the potential walls separating the different Chern-Simons numbers 
(and different vacua). The solution at the highest point of the potential is 
called a sphaleron, but also a transition that passes over this point is called 
a sphaleron (or a sphaleron process). These two processes will be outlined 
in the following sections .0 

^ Nonlinear coherent structures such as kinks, sohtons, sphalerons, vortices, and instan- 
tons arise in the study of a variety of interesting and important physical situations such 
as baryon number violation, properties of materials, optical communications, dynamics of 
phase transitions, and the physics of the early universe. 
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Electroweak 
potential 



Sphaleron process 
/ Sphaleron^ N 




a- Field 

1 Chern-Simons winding number 

Figure 1.2: Sketch of electroweak potential with minima, Chern-Simons 
number, sphaleron and instanton. 

Sphalerons 

Mathematically, a sphaleron is a special type of solution to a partial differ- 
ential equation. It is a saddle-point solution to the equations of motion with 
a single negative eigenvalue. 

It corresponds physically to the solution where the energy of the state 
is equal to that of the highest point of the (periodic) potential. For the 



electroweak potential, the energy of the spahleron can be calculated |12] to 
be 



Esph{T « 0) « 10 TeV. 



(1.7) 



The result is only an approximation, but it still gives an idea of the order 
of magnitude of the sphaleron energy. 

Note that the energy of the sphaleron actually depends on the temper- 
ature. The sphaleron energy decreases with increasing temperature, which 
means that when the temperature approaches the electroweak phase transi- 
tion temperature ~ 100 GeV, the sphaleron energy is believed to discontin- 
uously go to zero. 

The rate per volume of baryon number violating effects can be approxi- 
mated to be 



r(r) oc e 



-'sph 



{T)/T 
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for a certain regime |Q] (probably) below the electroweak phase transition. 
This is the same type of transition rate formula as for the pendulum, and we 
expect that near the critical temperature the exponential suppression factor 
will take effect [|l3| (when decreasing the temperature). 

As can be seen in the previous formula, the sphaleron processes will work 
against the production of a baryon-antibaryon asymmetry. In the scenario 
above, these sphaleron processes can be seen as the reason for the necessity of 
a thermodynamic non-equilibrium because the sphalerons vehicle the change 
in baryon number. 

Instantons 

Instantons are, like sphalerons, special types of solutions to partial differ- 
ential equations. Intuitively, they can be described as an energy configura- 
tion in space-time. Physically, they describe the tunneling between differ- 
ent Chern-Simons numbers/vacua. The tunneling rate betweeen different 
Chern-Simons numbers can be calculated to be 

Toce--^-^, (1.9) 

where Se is the Euclidian action obtained from the Minkowski-space action 
by a Wick rotation t — > —it = r. The tunneling rate at zero temperature 
can be calculated Q to be 

r(r = 0) ~ 10"^^°. (1.10) 

Thus, in our everyday life, instanton processes will play a very insignificant 
role, and baryon number is virually conserved. 

1.2.4 Charge and Charge-Parity Violation 

In the previous section a rather detailed study of the baryon number viola- 
tion has been done because it is closely related to the subject of this report, 
the hypercharge axion. As for the CP violation, it has nothing to do with 
the axion, and only some brief remarks will be made. 

Even though C and CP symmetries are known to be broken in the 
standard model, this source of CP violation is far too small to produce the 
observed asymmetry between matter and antimatter. 

In the minimal supersymmetric standard model (MSSM) CP violation 
can be enhanced by including a second Higgs doublet or by higher dimension 
operators, but this subject is too vast to be treated in this report. For a 
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resume of these effects, see This enhancement of the CP violation 
could be strong enough to produce the present dominance of matter over 
antimatter. 

1.3 The Axion 

The indroduction of the hypercharge axion has two objectives. First, it ex- 
plains the baryon number violation that gives rise to the matter-antimatter 
asymmetry in the universe. Second, it explains the existence of large-scale 
magnetic fields in diffuse astrophysical plasmas. The first Sakharov condi- 
tion is the principal subject of this study of the axion. The second Sakharov 
condition, the presence of C and CP violation is expected to be explained 
elsewhere and the third Sakharov condition is supposed to be fulfilled by a 
strong first-order electroweak phase transition.]^ 

The basic idea of the procedure from axion to baryon number violation 
is the following: The hypercharge axion couples to the hypermagnetic field, 
amplifying it to the extent that it gives sufficient baryon number violation 
to account for the matter-antimatter asymmetry. 

This section covers the theoretical background and properties of the 
hypercharge axion. 

1.3.1 Amplification of Hypermagnetic Fields 

In the universe today, we observe mysterious large-scale magnetic fields. 
These could originate from hypermagnetic fields that existed before the elec- 
troweak phase transition. The hypermagnetic fields are thought to behave 
basically like the ordinary electromagnetic fields, obeying Maxwell's equa- 
tions. With some different compositions there still exists two parts, one 
weak (SU(2)2,) and one hyperelectromagnetic (U(l)y). 

The part of the baryon number violating term, which is important for 
the hypercharge axion, is expressed as 

oc y^.y'^'' ocHy-Ey (1.11) 

because the hypercharge axion couples to the hyper(electro)magnetic fields, 
Hy and Ey- In the equation, Y^^i, is the hypermagnetic field strength and 
Yfif = ^fiuapy^^ is its dual {e^uap is the permutation symbol, defined in 
appendix 

^ According to long-range uniform magnetic fields could strengthen the electroweak 
phase transition to the point that the baryon number asymmetry could survive. 
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Normally, this source of baryon number violation is ignored because vac- 
uum is supposed in the initial as well as in the final statej^ However, in 
the theory of the hypercharge axion, we suppose that a hypermagnetic field 
exists before the electroweak phase transition. Where it comes from would 
be the issue of another study, but according to [l^ there "are neither com- 
pelling theoretical arguments nor motivated phenomenological constraints 
which could exclude the existence of magnetic fields prior to the nucleosyn- 
thesis epoch." The same article also gives the ratio for the baryon-to-entropy 
ratio 5{nB/s) due to the existence of the hyperelectromagnetic fields. The 
formula depends on space, time and hypermagnetic field strength. 

In [||] it is shown that the existence of a hypercharge axion introduces 
a new term in the Maxwell equations, which is a coupling between the 
hypercharge axion field X and the hypermagnetic field By- The solution of 
the Maxwell's equations can be decomposed in Fourier modes with different 
frequencies. The equation of motion for the hypercharge axion can be solved 
and the solution are periodic oscillations or rolling around the minimum of 
its potential. 

Depending on the frequency of the axion oscillation/rolling, it will cou- 
ple to (=come to resonance with) different Fourier modes. More specifically 
there will be a certain mode that will be maximally amplified, while the 
nearby modes are amplified somewhat less. The amplification is exponen- 
tial and can achieve values of 10^^ or larger only after a few cycles. Such 
high hypermagnetic fields can lead to enough baryon asymmetry genera- 
tion through equation ( |l.ll| ). In order to preserve the baryon asymmetry, 
the oscillation or rolling should take place just before or during the phase 
transition to avoid diffusion of the magnetic modes, which erases the field, 
and thus the conditions for baryon number violation. In the scenario of a 
rolling hypercharge axion, the amplification will affect very long wavelength 
modes. These modes do not diffuse as quickly as higher wavelength modes. 
This means that the resonance does not have to occur at the moment of the 
electroweak phase transition, but could occur somewhat earlier. 



1.3.2 The Lagrangian 

The Lagrangian density describing hyperelectromagnetic fields coupled to 
the heavy pseudoscalar axion X in the resistive approximation [16| of the 

Even when vacuum is supposed in the initial as weU as the final state the weak part 
will still give rise to a baryon number changing current [M. 
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highly conducting electroweak plasma is: 

L= (iv^xv^x - v{x) - iy^^y^. - J^y^ - ^xy^^y^,) 

(1.12) 

where Y^^, is the hypercharge field strength, V'^ is the hypercharge photon, 
g = det{g^u) is the determinant of the metric tensor, A ~ is a scaling 
constant, is the covariant derivative, V{X) is the axion potential and 
is the Ohmic current. The last term /x . . . represents the possibility that a 
fermionic chemical potential survives the unbroken phase of the electroweak 
plasma. Finally, e^^^ is the permutation symbol, which is 1 for all even 
permutations of e^^^, -1 for all odd permutations and for all other values 
of the indices. 

An explanation of the different terms in the Lagrangian can now be 
appropriate. 

^V^XV^X is the good old dynamic term, correponding to ^^^^ in classical 
mechanics. 

V{X) is the potential of the hypercharge axion, which can be approximated 
|16| to be ^MyX^, a harmonic oscillator potential. My is a mass-scale, 
see below. 

1 y/ii/y^^ is the hyperelectromagnetic energy, which is proportional to Ey ■ Hy 
[E and H being, respectively, the hyperelectric and hypermagnetic 
field). 

J^Y^ is the ordinary current of electric charge. 

— ^XY^^^Y^jj is the interaction between the hypercharge axion and the hy- 
perelectromagnetic field, see below. 

1.3.3 Couplings 

The part of the Lagrangian that is of interest for the detection in accelerators 
is the interaction between the hypercharge axion[| X and the hyperelectro- 

^ In fact, it is this coupling that has given the name to the hypercharge axion. The 
couphng is the same as for the quantum chromodynamics axion, which was proposed as a 
possible explanation for the non-existence of a strong equivalent to the weak CP violating 
term. For more information, see, e.g., |l4|. 
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magnetic field Yfj^^y: 

L = ^^V^-, (1-13) 

where Y'"' = e"^^^y„;3 is the dual of the hyperelectromagnetic field, My = 
^ is the mass-scale of the hypercharge axion, which in principle could be 
anything that is much larger than the mass of the axion, mx- We can 
obtain the Lagrangian in terms of Z^s and photons if we decompose the 
hypercharge fields into the Z^, and the photon, 

Y — 8 Y — 8 Y 

Y^ = Af, cos Ow - Z^sin9w, 

where is the photon field, Z^ is the Z^, and Oyy is the weak mixing angle. 
We obtain 

L = ^[ aid,A,d,A^) + b{d,Z,d,Z^) 

-c{d^A,dpZ,)-c{d^Z,dpA,)l ^ 

where a = cos^ Ow, b = sin^ Oy/ and c = sin 9y/ cos Ow- Here we can clearly 
see the interactions and their coupling constants. In the first term we see 
the X — > 77 coupling, in the second the X Z^Z^ coupling, and in the 
third and fourth the X — > Z^7 coupling. 



1.3.4 Branching Ratios 



The branching ratios of the processes can be obtained from the matrix ele- 
ments of the Lagrangian: 



1 



32'itM 



COS 



X~>ZOZO 



Y 
4 

327rM2 
1 



32itMy 



cos^ 6\Y sin^ Oyy 



(m2^-m|)3 



(1.16) 



where 9w is the weak mixing angle, mx is the mass of the axion. My is the 
mass-scale, and Mz is the Z^ mass. 

We note that the squared coupling constants from equation ( |1.15| ) are 
here as well as the mass-scale squared. The factor 4 for X Z^^ comes 
from the two last terms in equation ( |1.15 ) which are identical. 

The total width for mx ~ 1 TeV, My ~ 1 TeV is about 1 GeV. 
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1.3.5 Candidates for the Hypercharge Axion in Extensions 
of the Standard Model 



Brustein and Oaknin treat this subject in |17]. There are several pseu- 
doscalars in different extensions of the standard model that couple to the 
hypercharge as proposed in equation ( 1.13|) . In order to drive baryogenesis, 
a condition on the rolling/oscillating time must be fulfilled []T7| . 

One candidate would be the heavy Higgs pseudoscalar in the supersym- 
metric standard model, but this coupling is shown to vanish in the sym- 
metric phase of the electroweak theory. Other possible candidates, like the 
pseudoscalar component of sneutrinos, do not couple to the hypercharge. 

In string theory there are several possible pseudoscalars but the condition 
mentioned above turns out to be difficult to fulfil. Besides, the coupling to 
the hypercharge generally happens at temperatures much higher than the 
electroweak phase transition. In summary, it is possible that the hypercharge 
axion could be described by string theory, but the conditions are rather 
severe. 
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Chapter 2 



Phenomenology of the 
Hypercharge Axion 

This chapter discusses the experimental aspects and the phenomenology of 
the hypercharge axion. It describes the signatures and background for the 
detection of such a particle in the ATLAS detector at the future large hadron 
collider (LHC) at CERN. 

First the detector is presented briefly to situate the project in the con- 
text. Second, the characteristics of axion decay are discussed as well as the 
process that is the subject of this thesis 

q + q^q + q + -1 + -i ^q + q + X (2.1) 

with subsequent decay of the axion to photons and Z'^'s. Third, the approx- 
imations are presented. The interactions are simulated with a Monte Carlo 
program (PYTHIA) and the detector with ATLFAST. The semi-classical 
Weizsacker- Williams approximation is made for the emission of photons from 
the primary partons. Fourth, the general backgrounds (7, and jets) are 
discussed with their characteristics and their effect on our signals. Fifth, 
the different processes that originate from the existence of the axion are 
presented. There are several, but only a few can be detected in ATLAS. 

2.1 The Detector 

2.1.1 Lcirge Hadron Collider 

The large hadron collider (LHC) is a particle accelerator that is built at 
CERN (European laboratory for particle physics) and will be ready for ex- 
periment in 2005. The accelerator will bring protons and, later, ions into 



20 



head-on collisions with higher energies than ever before. These high energies 
win ahow physicists to probe deeper into the structure and the fundamental 
properties of matter and energy. The elevated energy will allow new types 
of particles to be created, detected and measured, but they will not be easy 
to detect. 

Along with the new, interesting physics there will be huge amount of 
background, or noise, coming from other processes that are already known. 
The principal reason behind this background is the fact that a proton is 
made of three quarks, which makes the collision a six-body problem. It is 
not easy to calculate and predict theoretically, nor to accurately detect all 
outgoing particles. 

LHC will be able to produce three basic types of reactions. First, it can 
collide proton beams, each proton having an energy of 7 TeV, for a total of 
14 TeV in the interaction. This will be the main type of collision at LHC and 
the luminosity]^ will reach 10^^ cm~^s~^, which is about 100 times higher 
than that of current colliders. Second, heavy ions, like lead, can be collided 
with a total energy of 1250 TeV in the collision, about 30 times higher 
than that of the Relativistic Heavy Ion Collider (RHIC) at the Brookhaven 
Laboratory in the US. Third, on long term it is planned to combine LEP 
and LHC. Electron-proton collisions can be obtained with total energies of 
~ 1.5 TeV, about five times higher than HERA at the DESY laboratory in 
Hamburg. 

2.1.2 The ATLAS detector 

The ATLAS (A Toroidal LHC Apparatus ) detector is one of the detectors 
on LHC. It is immense (like a five-story building) and is hermetically closed. 
The parts of ATLAS that are relevant for this research are: 

Inner tracker - measures the momentum of each charged particle. The 
inner tracker is inside a solenoidal magnet of 2 Tesla, bending the 
charged particle tracks and allowing measurement of their momenta. 

Calorimeter - measures the energies carried by the particles. There is 
one electromagnetic calorimeter and one hadronic calorimeter. The 
calorimeters are thick enough to let through only 10"'' percent of the 
entering particles. 

^ Luminosity £, is directly proportional to the number of particles in each bunch and 
to the bunch-crossing frequency, and inversely proportional to the area of the bunches at 
the collision point. 



21 



Muon spectrometer - identifies and measures muons. The muons in- 
teract so feebly witli matter tliat most of them pass the calorimeter 
without stopping and without causing a shower of particles. The muon 
spectrometer consists of chambers to measure their track. A strong 
toroidal magnetic field is present to measure their momenta. 

There is also a trigger system that selects hundreds of interesting events from 
thousands of millions of others and a data aquisition and filtering system to 
further classify and distinguish each event (collision every 25 ns). 

Obviously, the detector has a limited resolution, which, e.g., means that a 
quark will be observed as a shower, a jet. Also, photons with low transverse 
momentum will not be detected, and particles have to be separated by a 
certain angle to be identified as different particles. 



2.2 Characteristics of the Hypercharge Axion 
2.2.1 Interactions 

What is interesting for us is the interaction the axion would have with 
matter. Let us investigate the interaction of the pseudoscalar axion field X 
with the hypercharge field strength Y^^. This is described by the Lagrangian 

L = -^XYr.uY^"', (2.2) 

where My is a scale factor with the dimension of energy, and Y^^ = e'^^^'^Y^fs 
is the dual of the hypercharge field strength. This is the Lagrangian proposed 
by Brustein and Oaknin Obviously there are several other interactions 
in the system, including potentials and self- interactions, but these do not 
infiuence the results of this section. For a more detailed discussion, see sec- 



tion 1.3. We will take My as a free parameter and determine for which 
values of My, and for which masses of the axion mx, it will be possible to 
discover the axion in ATLAS . The actual values of My and mx depend on 
the theory invoked to make the connection between the matter- antimatter 



problem and the hypercharge axion, see section 1.3.5 . 

The hypercharge field-strength consists basically of Z^s and photons, 
which means that the interaction will be between the axion and either two 
photons, two Z^s or one Z^ and one photon. Hence, to identify the hy- 
percharge axion in the detector, we will look for these three types of decay 
products. With the particles coming out from the collision, we will look 
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for photons and^ Z^s and reconstruct the invariant mass^ of these particles. 
This mass should be equal to the mass of the axion, and with a lot of events 
it should be recognizable above the background. 

As the hypercharge axion is expected to be rather heavy, its decay prod- 
ucts will often have a high transverse momentum. Most other particles will 
more or less follow the direction of the beam due to the high momentum 
in the beam direction. Hence, another characteristic of the axion will be 
off-axis photons and Z^s. 

2.2.2 The Event 

Brustein and Oaknin proposed the channel f + f ^ Z*/'j* — > Z/'yX 
for detection of the hypercharge axion X. However, the cross section for 
this process is rather low and another possibility to detect the axion will be 
investigated here to complement their work. 

We will study emission of two photons (one from each proton) colliding 
to give the axion. In the proton, generally only one of the quarks (an up or 
a down quark) or a gluon will interact in the collision with another quark or 
gluon from the other proton. The event that is the subject of this research 
can therefore be written as (figure ^H] ) 

q + q^q + q + 'j + j^q + q + X, (2.3) 

where the axion can disintegrate as 

X ^7 + 7 

^ ^0 + 7 (2.4) 
^ + ZO. 

We can make a rough estimate of the number of events for a certain mass 
of the axion, say 1 TeV. The number of events is, by definition, the product 
of luminosity, time and cross section: 

Nev = Cxtxa 10^^ cm'^s-^ x 10^ s x 150 pb x 10"^^ b/cm^ 15. 

(2.5) 

^ In fact, we look for the decay products of the Z^s, see below. 
^ The invariant mass M is defined as 

M^=(E,i?.)'-iiE.p.ir, 

where Ei and pi are, respectively, the energy and the momentum of particle i, and — 
v'i +Vy + v1 is the squared norm of the vector v. 
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Figure 2.1: The event that is the object of this study. Two quarks within 
the proton emit two photons, which colhde and produce the axion, which 
disintegrates to photons and/or Z^s. 



Here the 77 ^ X cross-section, a, at the resonance is taken from with 
the mass-scale My = 1 TeV, and where the luminosity C is extrapolated 



from the graph in |18| and integrated over the resonance width of the axion 
(about 1 GeV). This is the luminosity for the production of two photons 
with invariant mass ~ 1 TeV during LHC operation at low luminosity. It 
is calculated assuming coherent emission from a proton. In fact, it can 
be much higher if the photons are emitted from quarks and the proton is 
allowed to disintegrate. The time t is taken to be one year of running, or 10^ 
s, and 10~^^ b/cm^ is a rough conversion factor from cm~^ to barn. This 
is detectable, though the signal will be rather feeble. This means that we 
have to be more careful with the background. To increase the signal we will 
assume that we run during three years of high luminosity at LHC, giving 
ten times as many events: 

Nev ~ 150. (2.6) 

The forward jets could in principle serve as a signature for the process, 
but the enormous jet background in the LHC will completely destroy this 
possibility. This jet background (also called QCD background) comes from 
interactions between quarks. These interactions produce outgoing quarks 
that will give showers of particles in the decay process. These showers are 
detected in ATLAS and are called jets. Because LHC works with proton 
collision (i.e., strong interactions) at extremely high energies there will also 
be a large number of jets, in the background. 

We can also use the fact that the decay products of the axion (the Z^s 
and photons) are produced virtually back-to-hack in the transverse plane, 
i.e., the angle between the two particles in the transverse plane is 180 de- 
grees. The reason for this is that the two photons producing the axion have 



24 



almost zero transverse momentum (in comparison with their longitudinal 
momentum) which comes from the fact that the quarks in the proton have 
very low transverse momenta. 

In summary, the characteristic of the axion are decay products X — > 77, 
Z^Z^ or Z'^'y with an invariant mass producing a bump, emitted back-to- 
back and with high -px- 

2.3 Approximations 

2.3.1 Monte Carlo Simulation 

The events have been generated in PYTHIA [|l9| which is a Monte-Carlo 
simulation program for particle interactions. To simulate the detector effects 
with its limits in resolution and acceptance ATLFAST has been used. This 
means, among other things, that photons with < 60 GeV will not be 
detected and that photons which are too close to the beam axis will be 
mistaken for jets. 

The generation of the matrix elements to calculate the interactions have 
been provided by Gilles Couture. He has used the Weizsacker- Williams ap- 
proximation for the emission of the two photons and the parton distribution 
functions from [^], mostly Appendix A. 

2.3.2 Weizsacker- Williams Approximation 

For the interaction between the two quarks, the Weizsacicer-Wiiiiams &p- 
proximation (also called equivalent photon approximation or almost-real- 
photon approximation) has been made. In principle the electromagnetic 
field of the charged particle is treated as (almost) real photons. These can 
interact with the almost real photons from the other charged particle. This 
approximation is only valid when the transverse momentum of the charged 
particle is virtually zero, and when they have ultra-relativistic speeds. This 
should be fairly well accomplished in our scheme, as the charged particle 
is a quark, and the quarks inside a proton have very little transverse mo- 
mentum and the speed is only a small fraction from that of light. Another 
consequence of this approximation is that the protons are supposed to split 
and produce fragments in the forward direction. 

There does not exist any good theoretical description of 77 processes 
from proton-proton collisions. However, the Weizsacker- Williams approxi- 
mation suffices well. 
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2.4 Background 



The background consists of all processes that give the same type of signal 
as the interesting event. An everyday example would be to weigh different 
animals, trying to distinguish them only by their weight. Suppose that we 
look for the elephant (the axion). We would have an enormous number 
of light animals, like insects, indistinguishable from each other, (loosely 
corresponding to the QCD background, see below). As the size goes up, 
fewer and fewer animals are heavy enough. Eventually, we might see a bump 
for humans around 70 kg, and around 5 tons we find the Indian elephant 
bump. However, there are other animals that could be this heavy, rhinoceros 
and whales for example, though it would be rare. This is the background. 

This section is dedicated to the different backgrounds for the axion. 
First, the interesting processes will be presented so that we know what 
backgrounds to consider. In summary, these backgrounds are processes pro- 
ducing photons and Z^s each of which will be discussed separately. A brief 
discussion of the jets is also included. 

In the calculation of the backgrounds for the different processes, the con- 
tribution of the Higgs boson is excluded. This is done because the Higgs 
boson has an unknown mass, which makes it difficult to estimate the back- 
ground it could produce for the process. In fact, we do not even know that 
it exists. If it exists as expected with a mass of 115 GeV, its contribution 
would be negligible, because the branching ratios of to photons and Z^s 
are very small [^]. 

2.4.1 Decay Channels 

The possible decay channels of the axion are (with branching ratios^ over 
the arrows): 

X > 7 + 7 

52% , „o 6% , r , 7 
> J + Z" > J + 1 + 1 



70% 



i% ^0 + ^0 036% i^l^i^l 



7 + jet + jet ^2 7) 



3.6% 
> 

36% 



1 + 1+ jet + jet 
jet + jet + jet + jet, 



4 



The branchi ng rat ios of the axion depend on its mass, but this dependence is very 



weak (see section 1.3.4). In the range 500 < mx < 5000 GeV the difference is less than a 



few tenth of percentage. 
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where we have used the fact that the decays to two specific leptons 
3 percent of the time, to two quarks (jets) 70 percent of the time and to 
neutrinos the rest of the time. The r is more difficuh to detect at the LHC 
and is not considered here, leaving the electron and muon, for a total of 6 
percent detectable lepton decays. 

As we can see, our primal branchings will be: 

X -^j + z^^j + i + i ^^■^> 

because X ^ + is rare and X ^ + jet + jet is very polluted by the 
QCD background. A more detailed discussion follows below in section 
Before discussing this further some general backgrounds will be described. 



2.4.2 7 Background 

We can have photon emission at different stages in the calculations. First, in 
initial and final state radiation. Second, in specific processes giving photons 
as end-products. Third, a jet can be taken to be a photon in the detector. 



Initial and Final State Radiation 

Initial and final state radiation are decays of particles that occur before 
and/or after the event. Examples of initial/final state radiation are g — > (?7 
and q qg. 



Photon Processes 

The principal processes that give rise to photon production in LHC are 



// - 




fg - 


- / + 7 


// - 


-> C/ + 7 


99 - 


g + j 


// - 


-^7 + 7 


99 - 


-^7 + 7 


// - 


- 7 + ^°, 







(2.9) 



where the fermion / is a quark in our case and g is the gluon. The processes 
f f ^ g + ^ and gg ^ g + ^ have a high cross-section, but on the other hand 
they do not correspond to the exact final state of any of our processes. As 
for the other processes, they all have a very small cross-section, leaving us 
with a fairly small photon background. 

Another process that could be expected to contribute is gg — > ^Z^ . The 
contribution would be to the background oi X ^ Z^j, but it is not included 
in PYTHIA and will not be included in our study. 
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Fake Photons 



The resolution of the photons is rather good, but in some rare cases (ap- 
proximately 1 time in 3000) a jet can be mistaken for a photon |23]. When 



the signal is strong, this is not a problem. However, when the signal is weak, 
this "fake" photon background could also be significant. 

2.4.3 Z° Background 

As two of the three branchings of the axion produce Z^s it is important 
to see which other processes give can be produced in several 

different ways: 

/ + / ^ 7/^° f + 9 ^ f + Z^ 

f + f ^ g + Z^ {Z^Z^ ^ ZO + ZO) 

/ + / ^ 7 + ZO {W+W- ^ ZO + ZO) ^^-^^^ 

/ + / ^ ZO + ZO 

Among these processes f + f g + Z^ and f + g ^ f + Z^ have high 
cross-section but do not directly lead to one of the final states of the axion 
decay either. This means that the missing photon will have to come from 
initital and/or final state radiation, as described in section p. 4. 2 . Fortu- 



nately, a high-pr photon from these sources is very improbable, diminishing 
the importance of these background processes. 

For Z°Z° ^ Z° + Z° and W+W' ^ Z^ + Z°, the cross-sections are 
extremely small, making it rather safe to exclude them from the calculations 
(see equation |37 



2.4.4 Jet Background 

The QCD processes that have been generated are 

qm 99 qi9 qi9 (2-ii) 

99 quqk 99 99, 

where the evolution and showers come from initial and final state radiation 
as described in section p. 4. 2 . 



In LHC there will always be an enormous background from jets due to 
the high energy and luminosity of the beam. This means that it will be very 
noisy to try to detect events with jets as final states. If the mass of the 
axion is relatively high, the criterion of high-pj- decay products can be used 
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to reduce the background significantly. Normally, the produced jets will 
follow the beam direction because the quarks have relatively low-px in the 
protons. However, the problem with a heavy axion is that its cross section 
for axion production decreases with increasing mass. This means that even 
with a low background, the cross-section for axion production might be too 
low to allow detection. 

We know that the signal studied should have two and only two jets. 
However, since we have used the Weizsacker- Williams approximation we 
cannot exploit this characteristic to apply cuts on the number and angle of 
the jets. 



2.5 Processes 



As can be seen from the branching ratios in the previous section and the 
estimate of number of events, equation (^^), there will be almost no events 
with production of the Z^s. Furthermore, LHC has a huge QCD background 
(jets) making it a difficult task to discern — > jet+jet. This leaves us with 
two principal processes, X^7 + 7andX^7 + Z°^7 + / + /. These wih 
be described in more detail below, while the others will be treated briefly in 
a separate subsection. 



2.5.1 Process X — i> 7 + 7 

This process should be very clean. Photons are rather easy to detect and 
their energies and momenta can be determined with good precision. Most of 
the noninteresting events (background) will be possible to sort out because 
the photons are produced back-to-back in the transverse plane, generally 
with high transverse momentum. 

The background processes that should be important are 

^ (2.12) 

gg -^11, 



along with the other processes in equation ( |2.9D with initial and final state 
radiation to produce the other photon. 

The 77 decay channel will be our principal source of events that hints at 
the existence of the axion because of its high branching ratio (44 percent) 
and the fact that photons are easily identified and (normally) do not decay. 
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2.5.2 Process 7 + ^° ^7 + ^ + ^ 

This process will also be very clean, but unfortunately also very rare. The 
total branching ratio for this process will be 3.1 percent and with the es- 
timated total number of events equation (|2.6|) this leaves about five events 
that will be produced during three years of high luminosity at the LHC. 
When we add the limits of the detector, along with the condition that three 
particles should be detected (7//), the number of observed events should be 
very small, at least for mx ~ 1000 GeV. 

The background processes that should be important are 



(2-13) 



and the other processes in equation ( p.lOD with initial and final state radia- 
tion providing the photon. 



2.5.3 Other processes 
Process X^7 + Z'^^7 + jet + jet 

The process X 7 + Z'' ^ + jet + jet has a fairly high branching 
ratio (36 percent), which gives a good number of such events. However, as 
mentioned before, the fact that the jets have to be used to reconstruct the 
invariant mass of the is very unfortunate. Even with the back-to-back 
criterion and the high-p-r criterion, the QCD background (see section 2.4.4) 
is overwhelming and effectively drowns the signal. 



Process X^Z^ + Z^-^l + l + l + l 

The process X Z^ + Z^ / + / + / + / should be very clean, but also very 
rare. Combining the low total branching ratio (0.0144 percent) with the 
difficulty for the detector to observe four distinct leptons leaves us without 
a single such event during three years of high luminosity at the LHC. 

Still, the event will be very clean because each pair of leptons can be 
combined to produce the invariant mass of Z^. The Z^s then give the 
invariant mass of the axion and, in addition, they should be back-to-back, 
like in the case of X ^ 77. 

Possible background events come from // — > Z^Z^. 
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Process X ^ + ^ I + I + jet + jet 

This process should be rather polluted because of the large QCD background 
that will give jets, and Z^ emission that will give the leptons, as discussed in 
section 2.4.3 . Furthermore, the signal in itself will be very weak (branching 
ratio 0.14 percent) and difficult to detect with its four particle final state. 
An example of background is 

t + t^ bW^ + bW'+ jet + jet + /+ + i/ + r + z^. (2.14) 



Process X ^ Z° + Z° ^ jet + jet + jet + jet 

The process X ^ Z^ + Z^ ^ jet + jet + jet + jet has only the back-to-back 
and high-pj' criteria to differentiate it from the enormous QCD background. 
These criteria are far too weak to make the signal discernable above the 
background. 
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Chapter 3 

Results and Analysis 



As discussed in the previous chapter, the simulation of the processes and 
the backgrounds has been done in PYTHIA with the axion interaction 
as an externally included process and with ATLFAST for simulation of the 
ATLAS detector. 

This chapter consists of four parts: The simulation of the signal of the 
different decay channels, the simulation of the background, both together, 
and finally some conclusions. The simulation of the signal is done for dif- 
ferent masses of the axion and gives us a hint of where a search could be 
fruitful. The signal will be stronger for lower masses, but the background 
will also be stronger. The analysis of the background will be done completely 
for mx ~ 1000 GeV. This gives an indication of which background processes 
are significant. This is important to know because at lower energies we have 
to generate more events as the cross-section goes up. The backgrounds that 
are found to be important are simulated also for mx ^ 800 GeV. 



3.1 Signal 

First of all we calculate the luminosity for different invariant masses of 77. 
The process is q + q ^ q + ^77. These are the photons that will interact 
to produce the axion and consequently their invariant mass should be the 
invariant mass of the axion. When we assume a total luminosity of 10^^ 
cm^s~^ we obtain the result shown in figure 3.1. Now, we calculate the 
cross-sections of the processes. With the cross-sections we calculate how 
many events we expect through 

TVe^ = £ X t X (3.1) 
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Figure 3.1: Luminosity of 77 pro- 
duction as a function of their in- 
variant mass for a pp luminosity 
of 10^'^ cm^^s^^. 



Figure 3.2: Number of events pro- 
duced for the different processes 
for three years of high-luminosity 
LHC operation. 



We will require at least ten events and we assume three years of high lu- 
minosity, which gives an integrated luminosity]^ of 300 fb~^. The number 
of events will be measured in units of (^^)^. However, we lose more than 
half of the events in the ATLAS detector due to problems with resolution 
and the fact that the detection of photons is not possible at near forward 
angles. This means that we need at least A'^e,; > 20 to really detect the 
particle. Figure |3.2| shows the number of axions produced as a function of 
axion mass, mx, for the different decay channels. We see that for an axion 
mass mx ~ 1000 GeV the number of events is almost a factor ten higher 
than the one predicted from equation (|2.6| ). The reason for this is that the 
photon emission is done directly from the quarks, and the proton is allowed 
to disintegrate (see explanation following equation ( ^.51) ). Note that this is 
only a first estimate and that the background has to be considered too. This 
is done in the following section. 



^ Integrated luminosity is £ x t. 
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3.2 Backgrounds 



In this section we will discuss preselection cuts on the background, which 
means that we limit the kinematic phase space in order to decrease the 
cross-section. Otherwise we would have to generate an enormous amount 
of background. We also generate the complete backgrounds for mx ~ 1000 
GeV and the significant backgrounds for mx ~ 800 GeV. The signal and 



the backgrounds are discussed together in the next section, 3.3 



3.2.1 Cuts 

In order not to a have to generate an overwhelming number of background 
events we introduce preselection in the calculations. A preselection means 
that an event is rejected at an early state because it does not fulfil cer- 
tain criteria. The cuts in our simulation have been made on the transverse 
momentum px and the center-of-mass energy of the two partons par- 
ticipating in the reaction. As discussed in chapter || a characteristic of the 
axion is high-p-r decay products, motivating a cut exclude pT < 300 GeV 
which should be valid for mx > 700 GeV. In order to even create something 
similar to the axion, the center-of-mass energy should be comparable to the 
mass of the axion. We have chosen to exclude \/I < mx — 100 GeV. 

To validate these cuts, the backgrounds have been simulated with lower 
cuts. With proper normalization^ the two signals are then compared. In the 
region around mx the shape of the curves should be similar. An example is 



shown in figure 3.3 



These were preselection cuts, which act in the generation of events to 
reduce the cross-sections. In consequence of the 3^^^ characteristic of the 
axion we apply a third cut, but this can only be done after the generation 
of events. This is a cut on the angle in the transverse plane between the 
particles that could have been the decay products of the axion. How these 
particles are selected is discussed below. The cut on the angle is 3.0 < cp < n. 

Obviously, all these cuts must also be applied on the signals themselves. 
This will reduce the signal, but even more so the background, which in 
principle should be randomly distributed. 



Selection Criteria 

The selectrion criteria for the different processes are: 

^ When two signals have different cross-sections they must be normahzed. The nor- 
maUzation factor is = f^T^' where Ni is the number of generated events for signal i. 
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Figure 3.3: Comparison of reconstructed m^^ from background ff^j + ^, 
with preselection cuts, \/i > 900 (full histogram) and V§ > 700 (dashed 
histogram) . 



X — ^ 77 Select the two photons with highest transverse momentum and 
calculate their invariant mass. 

X — Z°7 — 11^ Select the two leptons with highest transverse momentum, 
calculate their invariant mass (should give the Z'^), select the photon 
with highest transverse momentum and calculate the invariant mass 
with the reconstructed Z^. 

X — > Z^j jet + jet + 7 Same as Z^ II except that the leptons are 
replaced by the jets. 

X — > Z^Z^ nil Select the four leptons with highest transverse momen- 
tum, calculate the invariant masses between all the leptons, choose 
the pair that gives the mass closest to that of the The other pair 
is supposed to give the other Z^, and together the pairs give the mass 
of the " axion" . 

X Z^Z^ II + jet + jet Same principle as for llll. 

X Z°Z° ^ jet + jet + jet + jet Same principle as for Z°Z^ llll. 
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QCD background In some rare cases (about 1 in 3000) a jet can be mis- 
taken for a photon and the selection process for this is to choose the 
photon and the jet with highest pT and calculate their invariant mass. 



3.2.2 Complete Background for nix ~ 1 TeV 

As stated in the previous chapter, there are several different backgrounds 
that will affect the signal. The cross-sections are calculated with the two 
initial cuts, px > 300 GeV and \/I > 900 GeV. For three years of high 
luminosity at the LHC, the integrated luminosity is 300 fb~^ and the number 
of expected events is calculated from this. The specific backgrounds that 
have been treated in more detail are: 



Initial 


Final 


Actual # of events 


Simulated # of events 


// 




46200 


59910 


// 


g + 1 


202020 


400000 


// 


- g + z^ 


303300 


400000 


// 


-^7 + 7 
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10000 


// 
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// 


ZO + ZO 


6732 


34410 


fg 


- / + 7 


886500 


400000 


fg 


- f + z^ 


1261500 


400000 


z^z^ - 


z'^ + z'^ 








WW 


z'^ + z^ 


3 





gg 


-^7 + 7 


66 


10000 


gg 


g + 1 


57 


10000 


7 + jet - 




564450 


400000 








(3.2 



Here actual 7^ of events is the number of events that would be produced for 
an integrated luminosity of 300 fb~^. Note that the last line represents the 



fake photons that come from QCD processes, as described in 2.4.2 . As can 
be seen the Z^Z^ Z^ + Z^ and WW Z^ + Z^ have a very small cross- 
section. In fact, they are so tiny for the given cuts (see section 3.2.l|) that 
statistically, there will not be a single such event detected in ATLAS. These 
events are excluded as they will not contribute to any of the backgrounds. 
For most processes, an excess of events has been produced. The exceptions 
are fg—>-f + l and fg—^f + Z^ and the photon-jet confusion in the QCD 
background. However, this mixup of photons and jets is only considered 
for the X — > 77 channel. That these processes are underproduced can be 
motivated by the small contribution they add to the total background. 
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Figure 3.4: Histograms for the contributing backgrounds for the process 
X — > 77 for nix ~ 1000 GeV, after preselection cuts. 



For the sake of completeness, all these backgrounds have been treated 
with all the selection criteria of the different processes of interest, see section 



13. However, only a few of them will actually produce a background for a 



certain signal. These are the backgrounds presented in the figures. 
Background for X ^ 77 at nix ~ 1 TeV 



The background of this process is presented in figure 3.4. The signals are 
normalized to an integrated luminosity of 300 fb~^. We notice that the only 
processes that actually give a background are // ^ 7 + 7 and 55 — > 7 + 7- 
Even though their cross-sections are not so high, they have a final state 
that is irreducible. Hence, these are the two processes that we will keep 
when we calculate the backgrounds for different axion masses. Of these two, 
// ^ 7 + 7 is largest, which is only to expect in proton-proton collisions. 

In addition to the direct backgrounds we have also simulated the QCD 
background to account for the possibility of confusion of a jet and a pho- 



ton, as discussed in chapter p. 4. 2 . This is simulated just like all the other 



processes except that the cross-section is taken to be 3000 times^ less (other- 
This is motivated by the result, which does not contain much background (figure I 
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Figure 3.5: Histograms for the contributing backgrounds for the process 
X Z^^ — > /Z7 for mx ~ 1000 GeV, after preselection cuts. 



wise the number of events would be enormous). When the selection criteria 
are apphed, the invariant mass is calculated from the photon and jet with 
highest transverse momenta. 



Background for X —^ jZ^ — > 7// at mx ~ 1 TeV 



The background of this process is presented in figure |3.5| . The signals are 
normahzed to an integrated luminosity of 300 fb~^ and the cuts (see section 
3.2.1 ) have been applied. The only backgrounds that contribute to the signal 
are // — > ^/Z^ and // 'y + Z^. These two signals will be kept for analysis 
of different axion masses. We do not consider the possible confusion of a 
jet with a photon because the probability of production of a Z^ with the 
imposed cuts is very low in these processes. In analogy with the X — > 77 
process, the // — > 7 + Z'^ gives the irreducible background. 



Background for Other Processes at mx ~ 1 TeV 

As predicted in chapter |2|, the background for axion decay to jets is tremen- 
dous, several order of magnitudes larger than the signals themselves. The 
one with lowest background (but still very large) is X ^ Z°Z° jet + jet + 
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Figure 3.6: Histograms for contributing backgrounds for X ^ 77 at mx ~ 
800 GeV, after preselection cuts. 



//, but unfortunately the cross-section is too small. Thus all final states that 
include jets are worthless for detecting the axion. The only process left, 
X Z^Z^ nil has a vanishing cross-section and will give about one 
event in 55 years of LHC operation at high luminosity (total: 6000 fb~^). 



3.2.3 Significant Background Processes for nix ~ 800 GeV 

The cross-sections are calculated with the two initial cuts, pT > 300 GeV 
and VS > 700 GeV (see section 3.2.1). The signals are normalized to an 
integrated luminosity of 300 fb~^. 



Background for X — > 77 at mx ~ 800 GeV 



As we found in section 3.2.2 there are only two important processes for 77 
production. However, we have also considered the QCD background as this 
is known to increase dramatically for lower cuts. The result is presented in 
figure 3^. Except for the lower cutoff in it is the same background as 
for the case of nix ~ 1000 GeV. 
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Figure 3.7: Histograms for contributing backgrounds for X Z^^ — > 11^ at 
mx ~ 800 GeV, after preselection cuts. 

Background for X 7^° -/II at mx ~ 800 GeV 



As we found in section |3.2.2| there are only two important processes for 7// 



production. The result is presented in figure |3.7| . No spectacular deviations 
from the case mx ~ 1000 GeV is noted. 

Background for Other Processes at mx ~ 800 GeV 

For completeness, the other possible final states for axion decay should be 
considered also for this cut. Due to the weak signal of Z^Z^ production and 
the enormous background for X ^Z^ j+jet+jet the above mentioned 
processes are the only ones that can be used for detection and identification 
of the hypercharge axion. If the integrated luminosity and the energy could 
be drastically increased they could be interesting, but in LHC they would 



hardly show up at all, as discussed in section 2.5.3 



3.3 Detection of the Hypercharge Axion in ATLAS 

In this section we will study the two main decay channels of the axion, 
X ^ 77 and X 'yZ^ — > jll. We will study the signals with the cuts 
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and apply the backgrounds to see whether the axion can be expected to be 
discovered. The condition for discovery will be that the signal should contain 
more events than five times the standard deviation of the background (i.e., 
five times the square-root of the number of background events in the mass- 
window), and that the number of events from the signal is > 10: 

Ne^>5xVB and N^^ > 10. (3.3) 

The ratio Nev/V^ is called significance. The number of events is integrated 
between mx — 20 GeV and mx + 20 GeV. Note that the signal is expressed 
in units of (^^)^, which means that the signal will go down rapidly with 
increasing My. Note that often in the theoretical models (see chapter |l|) 
My > Mx. 

Two different masses of the axion are discussed, 800 GeV and 1000 GeV. 



3.3.1 Signal and Background for X ^ 77 

For both masses, the signal is rather clean, and pretty strong. 

For mx = 1 TeV, the signal and the backgrounds are presented in figure 



3.8 and we can clearly see the signal over the background. In the window 
980 < mx < 1020 GeV we have about 204 background events and 210 
signal events. This gives a significance of 14.7 sigma with My = 1 TeV (see 
beginning of section). 

For mx = 800 GeV the number of events is higher but so is the back- 



ground, as can be seen in figure 3.£. In the window 780 < mx < 820 GeV we 
have about 405 background events and 255 signal events. The significance 
is about 12.7 sigma, which means that the signal is possible to detect with 
My = 1 TeV (see beginning of section). We remark that the background 
grows faster than the signal as the mass goes down. 



3.3.2 Signal and Background for X 7^° 'jll 

This signal is also clean. The only drawback is that the branching ratio for 
II is so low. 

The background is very small for mx = 1 TeV but so is the signal as can 
be seen in figure ^l0| . In the window 980 < mx < 1020 GeV we have about 



26 background events and 11.4 signal events. The signal is hardly above the 
background and we have a significance of 2.2 sigma for My = 1 TeV (see 
beginning of section). 

For a lower axion mass, mx = 800 GeV, the signal is not clearly visible 



either (see figure 3.11 ). In the window 780 < mx < 820 GeV we have about 
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Figure 3.8: Histograms for signal and background for X 77 at My = 1000 
GeV, mx = 1000 GeV, after preselection cuts. 
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Figure 3.9: Histograms for signal and background for X ^ 77 at My = 1000 
GeV, mx = 800 GeV, after preselection cuts. 
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Figure 3.10: Histograms for signal and background for X — ^ Z^j — > /Z7 at 
My = 1000 GeV, mx = 1000 GeV, after preselection cuts. 

45 background events and 15 signal events. The significance is 2.2 sigma for 
My = 1 TeV (see beginning of section). 

3.4 Discussion and Conclusions 

In the previous sections we have seen that there are basically two processes 
that can be used for detection of the hypercharge axion; X — 77 and 
X Z^^i llj. The former signal has the higher signal-to-background 
ratio. In the case of mx = 800 GeV and My = 1000 GeV the axion would 
be observed (a 14.7 sigma significance). The other signal, X Z^'y 
11^, is more difficult to distinguish and for the same parameters we have a 
significance of 2.2 sigma. This is not enough for discovery, in the case of a 
clear signal for the X — > 77, but it could support a discovery and verify the 
branching ratios. 

If we want to be certain that it is really the hypercharge axion that is 
detected, the branching ratios of X ^ 77 and X Z^'y 11^ have to 
be compared. The ratio T^^/Vii^ should be 28 for the hypercharge axion. 
The cross-section in itself cannot serve as a certain identification, though 
well as an indication, because the cross-section is proportional to l/My and 
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Figure 3.11: Histograms for signal and background for X Z^^ llj at 
My = 1000 GeV, nix = 800 GeV, after preselection cuts. 



My is an unknown parameter with dimension of mass. If all three decay 
channels are discovered, then it is a strong indication that we have discov- 
ered a pseudoscalar with the right coupling to the hypercharge. At that 
stage, the mass-scale My can be used to determine theoretically whether 
the pseudoscalar could amplify the hypermagnetic fields or not |17|. 

The scenario of mx = 800 GeV and My = 1000 GeV is taken from 
the article by Brustein and Oaknin |^], which gave the original idea of this 
thesis. The cross-sections in both their and our case are about 10 fb~^. 

There is also a theoretical argument why a particle in the mass region 
300 < mx < 1000 GeV should be the hypercharge axion. It is not expected 
that there exist any particles in this mass range with electroweak couplings. 
The Higgs could in principle have this mass but this is outruled by the LEP 
(Large Electron Positron collider) high precision electroweak measurement. 
Furthermore, for a Higgs of this mass, the decay channel to two photons 
would be non-existent anyway. The Higgs would also have other, stronger, 
decay channels (like 

Note that in the calculations in this chapter. My is supposed to be 1 
TeV, or, in other words, everything is expressed in terms of (^^)^- The 
parameter My is an energy scale, which is supposed to be much larger than 
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mx but in this thesis we disregard this to discover the hmits of hypercharge 
axion detection in ATLAS. For My < 1.6 TeV detection should be possible 
but hardly much above that. The precise limit depends on the mass of the 
axion too. 

This thesis covers the possibility of discovering the hypercharge axion in 
ATLAS, but there arc several points that could be elaborated and treated 
in more detail. Examples are the Weizsacker- Williams approximation, the 
parton distribution functions and the amount of background produced. A 
study of processes and backgrounds could also be made for other values of 
the hypercharge axion mass. 

There are also several questions in the theory behind the hypercharge 
axion that are still obscure. For example where the original hypermagnetic 
field came from and how enough CP violation can be obtained. 
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Appendix A 



Glossary 



Anomaly An anomaly is the failure of a symmetry to survive renormaliza- 
tion. 

Antimatter Antimatter is constituted of antiparticles. 

Antiparticle An antiparticle is defined as having the opposite quantum 
numbers as the corresponding particle, but the same mass. Parti- 
cles with their quantum numbers zero, except for spin which should 
be integer, (like the photon) are their own antiparticles. Experimen- 
tally there is no evidence for the effect of gravity on antimatter, but 
according to Einstein's principle of equivalence there should be no dif- 
ference to matter. However, there are other possible theories, cousins 
of Einsteins theory of gravitation, that predict weaker gravity, or even 
inverse gravity for antiparticles. For a review of antiparticles, see, e.g., 
J. Eades and F. J. Hartmann, Rev. Mod. Phys. 1999 71, 373-420. 

Axial current The axial current is defined as 

= ^/;(x)7^75V'(x), (A.l) 

where ^p{x) is the field at space-time point x, e(x) is an arbitrary 
function of x and 7q, (a = 1, 2, 3, 4, 5) are defined in appendix ^ 

Baryogenesis The generation of baryon asymmetry. 

Branching ratio The relative probability that a particle will decay in a 
specific way. Example: — > e^e~ has a branching ratio of about 3 
percent. The sum of all branching ratios should be unity. 

Charge (C) C is an operation that reverses the charge of a particle. 
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Chern-Simons number A wiiiding'iiuinber for the electroweak Lagrangian. 
The change in Chern-Simons number is proportional to the change in 
baryon number, (see section 1.2.1). 



Decay channels The possible decays of a particle. 
Example: 

— > nn, (id, ss, cc, 66, tt, e~e^, r^r^, feZ^e, Vt^ti 

where the different processes have different branching ratios. 

Electroweak phase transition The electroweak phase transition is when 
the electroweak symmetry is broken and the Higgs particle obtains 
a vacuum expectation value. It is estimated to have occurred about 
10"^" s after big bang, when the temperature of the universe was ~ 100 
GeV. 

Final state radiation Final state radiation is decays of particles or emis- 
sion of gluons or photons that occur after the event. In PYTHIA, this 
is the process that allow particle showers to evolve. 

Goldstone boson A zero spin, zero mass particle associated with symme- 
try breaking. 

Higgs boson A hypothetical particle required in the standard model of 
particle physics. The Higgs boson explains why W^, and have a 
mass. 

Initial state radiation Initial state radiation is decays of particles or emis- 
sion of gluons or photons that occur before the event. 

Massshell On the mass shell means that the particle is real (M^ = E'^ — 
IpI^). Off the mass sheh means that it is virtual (M^ / — 

Parity (P) Parity is an operation that gives a particle the same properties 
as if is was observed in a "point-like mirror" . In other words, the spin of 
the particle will be inversed. P {parity = p) = {parity = —p). The left- 
parity operator is represented by Pl = ^(1 — 75), where 75 is defined 
in appendix^ The right-parity operator is defined as Pr = |(1 -I-75). 

Sphaleron An unstable solution of a partial differential equation. A sphaleron 
process is the process of passing a spahaleron, (see |1.2.1| .) 
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Symmetry A symmetry operation does not change the physical solution. 
For example, the position of the origin in the coordinate system does 
not change the physical solution to problem. 

Winding number A winding number is a characteristic that does not 
change during continuous transformations, if a certain point is avoided. 
A trivial example of this is found in complex analysis: 



where the value of the integral is independent of the curve jab ^ long 
as the points A and B are the same and p is avoided. 

The winding number is a sort of shortcut which lets us know the 
consequences of micro-level rules without elaborate calculations. 

Quantum numbers The numbers that can be said to best describe the 
state of a particle. Examples: electric charge (Q), lepton number (L), 
baryon number (B), parity (P), spin (S), isospin (7), strangeness (S), 
and charge conjugation (C). 




(A.2) 
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Appendix B 

List of Symbols 



B.l In Equations 



X 


Pseudoscalar axion 


mx 


Mass of the axion 


My 


Mass scale 


Y 


Hyperchargc field strength (from U{1)) 


Y 


The dual of Y^,, = e^p^.Y^^ 




The permutation symbol, defined as +1 for even permutationes of e' 




-1 for odd permutationes, and for other values of the indices 


A, 


Hyperphoton vector 




vector 


mz 


Mass of 


H 


Hamiltonian 


T 


Temperature 


E 


Energy 


B = nB 


nB = — n^, baryon number: 




difference between the number of baryons and antibaryons 


L 


Lcpton number 


s 


Entropy 


r 


Transition rate or branching ratio 


Ncs 


Chern-Simons (winding) number 


Se 


Euclidian action 


Hy 


Hypermagnetic field 


Ey 


Hyperelectric field 


Jb 


Baryon current 




Chemical potential 
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7 
I 

9 
7 
/ 
Q 

jet 
t 
b 
P 

PT 

L 
a 
t 



Photon 
Lepton 
Gluon 
Photon 

Fermion (quark or lepton) 
Quark 

Jet of particles, produced by a quark or a gluon 
Top quark 
Bottom quark 
Momentum 

Transverse momentum, pT = 

Luminosity 

Cross-section 

Time 

Number of families of quarks/leptons, believed to be three 
Imaginary number, defined as i' 



Px ~^ Py if beam direction 



7 matrices 70 



71 



r 1 1 




010 


75 = 


00-10 


.00 -1. 




- 1- 




010 


72 = 


0-100 


. -1 00. 





10 
1 
10 
10 

i 
i 
-j 



= -1 






r 01 on 


73 = 


000-1 


-1000 




L 1 J 



B.2 List of Constants 

9w 28.7 degrees The weak mixing angle sinOw ~ 0.23 

h 1.054571597xl0"^'* Js h = h/lir, where h is Planck's constant 

c 299792458 ms~^ Speed of light in vacuum^ 

ks 1.3806503x10^23 jj^-i Boltzmann's constant^ 

eV 1.6022x10'^^ J Electron Volt 



Unless otherwise specified, we set h — c — ks = 1 
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B.3 Abbreviations 



ATLAS 


A Toroidal LHC Apparatus 


ATLFAST 


A FORTRAN program to approximately simulate the ATLAS-detector 


BAU 


Baryon Asymmetry of the Universe 


C 


Charge-reverse operator 


CERN 


European laboratory for particle physics 




(earlier: Conseil Europeen pour la Recherche Nucleaire) 


CP 


Charge-Parity operator 


CPT 


Charge-Parity-Time operator 


EW 


Electro Weak 


EWPT 


Electro Weak Phase Transition 


HCA 


HypcrChargc Axion 


LHC 


Large Hadron Collider 


LEP 


Large Electron Positron collider 


LTU 


Lulea University of Technology (in Swedish: Lulea Tekniska Universitet) 


MSSM 


Minimal Supersymmetric Standard Model 


P 


Parity operator 


PYTHIA 


A FORTRAN program to simulate collisions between particles 


QCD 


Quantum ChromoDynamics, the field theory for strong interactions 


QED 


Quantum ElcctroDynamics, the field theory for electromagnetic interactions 


SUSY 


SUper SYmmetry 


T 


Time operator 
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